Formulation of thermodynamics for the glassy state: 
configurational energy as a modest source of energy 
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Glass is an under-cooled liquid that very slowly relaxes towards the equilibrium crystalline state. 
Its energy balance is ill understood, since it is widely believed that the glassy state cannot be 
described thermodynamically. However, the classical paradoxes involving the Ehrenfest relations 
and Prigogine-Defay ratio can be explained when the 'effective' or 'fictive' temperature of the slow 
modes is taken as additional system parameter. 

Having straigtened out the proper picture, it is interesting to reconsider glass from a thermody- 
namic viewpoint. A shard of glass, kept at fixed temperature and volume, relaxes towards lower 
energy. Heat is released, inducing apparent violations of all basic thermodynamic laws. The most 
interesting application is to use glassy or amorphous systems as a source of energy, by extracting 
the configurational energy in a process of melting followed by crystallization. 
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Glasses have been of use to mankind from early on, 
be it as arrowheads for the stone age people of Corsica 
and the Americas, the obsidian battle axes and swords 
of the Aztecs, pumice scrappers for animal hides, or the 
tektite ornaments and fertility symbols of our ancestors. 
It has also been a source of fascination once its manufac- 
turing process became understood around 2000BC in the 
Caucasus and has continued to amaze us when consider- 
ing artifacts from the tomb of Thutmose 111(1480 BC), 
stained glasses windows in Chartres or the exciting world 
of sihca glass based fiber optics. 

Besides silica based ones, there exist glasses such as 
borate, phosphate, or germanate glasses, which have as 
primary and sufficient attribute the four fold oxygen co- 
ordination of the glass forming cation. Glassy ice can 
exist as crystals in forms topologically equivalent to sil- 
ica polymorphs. Also oxygen free glasses as the chalco- 
genide or halide ones exists, as do polymer glasses which 
go under the house hold names of nylon, PVC bottles, or 
handy wrap. For more examples, see Q. 

The archetypical glass for most of us is window glass, 
a fairly late invention around 500 AD, which commonly 
consists for about 70 percent of silica, the remainder be- 
ing partitioned in about equal proportions between soda 
and lime. It forms films which flow slowly causing a thin- 
ning at the upper end in the course of time, and time 
lapse pictures taken every 500 years would reveal that its 
behavior would be akin to that of a soap film. 

What all these glassy systems have in common is that 
they will structurally relax from an excited glassy state 
to a crystalline ground state in due course. Such relax- 
ation may be slow as is typical for systems with broken 
ergodicity, i.e. systems which do not reach equilibrium 
within a characteristic observation time. 

It is the nature of this structural relaxation which we 
wish to consider here in the light of a recently developed 
thermodynamic structure for glassy systems. After doing 
that, we shall reconsider the glass from a thermodynamic 



view point, and find apparant violations of all basic laws. 



The glassy state: a diverging time scale 

Any liquid cooled sufficiently fast will become glassy, 
i.e. it will lack the time to evolve into a long-range 
ordered crystalline array. In such glassy materials two 
types of processes can be distinguished: fast or fi pro- 
cesses, which happen on short scale and remain in ther- 
mal equilibrium, and slow or a-processes, which are 
restructurization processes that involve many atoms. 
Somewhere during cooling they fall out of thermal equi- 
librium, or, more precisely, their relaxation time exceeds 
the time needed to reach the (sliding) equilibrium. Re- 
cently the microscopic nature of the liquid landscape that 
leads to these slow processes was discussed ^ 0. 

The relaxation time in glassy systems scales with the 
viscosity, which is a measure of the internal friction of a 
system. It has in general a strong temperature depen- 
dence varying in silica-based glasses commonly over 18 
orders of magnitude in Pa.s, and is the least understood 
of all glass properties. Both viscosity and relaxation time 
will grow sufficiently rapid in a typical cooling experi- 
ment that a point is reached where equilibrium no longer 
exists between the thermal state of the system and the 
surrounding temperature. This point is called the glass 
transition temperature, commonly occuring in the tem- 
perature range of roughly two thirds of the melting tem- 
perature in silica based glasses. 

The growth in relaxation time for so-called frag- 
ile glasses, such as toluene, glycerol, and ethanol, is 
built into the Vogel-Fulcher-Tammann-Hesse law Tcq ~ 
exp[S/(r - To)] § H §. Adam and Gibbs § pro- 
vided an intuitive explanation of this law and Parisi |^] 
gave a field theoretical argument for the related behav- 
ior Tcq ^ exp{[B/(r - Tq)Y~^}, where d is the spa- 
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tial dimension of the system (d = 3 in practice). For 
so-called strong glass, such as Si02 and Ge02, on the 
other hand, the relaxation time follows an Arhenius law 
Tq, exp{A/T). Somewhat related expressions were pro- 
posed by Kivelson et al. ||] and by Schultz |0 to 
describe relaxation times for all known glasses todate. 

In our studies we have benefited from exactly solving 
simple model systems with an Arrhenius law 
lUl. or a VFTH-law (l|Jl|l. In models such as soft-spin 
models the spherical p-spin model jT^, or a 

directed polymer model |19[| partial answers can be ob- 
tained. Together this has resulted in a coherent picture 
that we discuss in the following. 



Breakdown of thermodynamics for the glassy state? 

Attempts for a thermodynamic interpretation of a 
glass phase go back to the twenties. Let us consider the 
experimental way to determine the entropy: one inte- 
grates d^* = {Cp/T)dT, where Cp is the measured heat 
capacity, from a sufficiently high reference temperature in 
the liquid phase, down to lower temperatures. The excess 
entropy is the thus obtained entropy minus the entropy 
of the crystalline state. It decreases with temperature, 
but will remain positive down to zero temperature, thus 
resulting in a residual entropy, in apparent violation with 
the third law. 

The glass transition temperature depends on the cool- 
ing rate. Slower cooling provides the system with more 
time to relax, with as result that it will remain in equi- 
librium to lower temperatures. Kauzmann made it plau- 
sible that, for infinitely slow cooling, AS* goes to zero 
near a specific temperature, now called the Kauzmann 
temperature Tk |2|]. It is to be expected that Tk and 
To coincide, and this indeed occurs in an exactly solv- 
able model glass ||lj,|l5| . The ultimate vanishing of the 
excess entropy near this point is called the Kauzmann 
paradox. A thermodynamic Kauzmann transition has 
been found in mean field models, like the Random En- 
ergy Model the spherical p-spin model and 
a polymer model p^ , but also in a model obeying a 
VFTH-law |l^Jl^. Finally, it also shows up in a mean 
field approach (so-called replica approach) to the glass 
transition in realistic models for glass [ pO[ . 

Proposals to formulate properties of the glassy state in 
a thermodynamic frame have remained unsuccessful till 
recently when we adressed this issue. The negative con- 
clusions were not based on subtle entropy arguments, but 
on testing the so-called Ehrenfest relations (see the Ap- 
pendix) . They are relations between the smeared out dis- 
continuities that occur around the glass formation tem- 
perature in cooling experiments at constant pressure, in 
the specific heat Cp — d{U+pV)/dT\p and the expansiv- 
ity a — d\nV/dT\p. One can treat them as true discon- 
tinuities, because of the small size of the smearing region. 
When comparing cooling experiments at different pres- 



sures one also finds discontinuities in the compresibility 
K — —dlnV/dplr- A glass transition seems therefore to 
be a textbook example of a mean field second order phase 
transition. However, this analogy is not perfect because 
there is no diverging length scale, and the smaller value 
of the specific heat occurs in the low temperature regime. 

The two Ehrenfest relations are particularly useful 
tests for glasses. The common notion of half a cen- 
tury of research is that the first relation is always vi- 
olated, whereas the second relation is obeyed in most 
cases, though not in all |2^. In the most careful 
experiment we came across, the glass transition in atac- 
tic polystyrene with a cooling rate of 18 K/h at ambient 
pressure, ACp is too small by a factor 0.77 |2^], while 
H = 1.09 « 1. 

The failure of in particular the first Ehrenfest rela- 
tions has led to the general belief that "thermodynamics 
does not work for glasses, as there is no equilibrium". 
This unapplicability to glasses is unsatisfactory on two 
accounts. Firstly, thermodynamics should also apply to 
non-equilibrium systems. Secondly, in view of the large 
time scales, up to many years, one would hope for a uni- 
versal description of the energy balance. 

The negative conclusion about the applicability of 
thermodynamics was mainly based on the failure to un- 
derstand the Ehrenfest relations and the Prigogine-Defay 
ratio. It should be kept in mind that, so far, the ap- 
proaches leaned very much on equilibrium ideas. Well 
known examples are the 1953 Davies- Jones theory [ p7[ , 
the 1958 Gibbs-DiMarzio H and the 1965 Adam- 
Gibbs § papers, while a 1981 paper by DiMarzio has ti- 
tle "Equilibrium theory of glasses" and a subsection "An 
equilibrium theory of glasses is absolutely necessary" . 
In particular, the Gibbs-DiMarzio model has long been 
considered as an important theory for glass. It is our pur- 
pose to explain that the step to an equilibrium theory has 
been an unfortunate move. 



Thermodynamics far from equilibrium and the 
effective or fictive temperature 

Though a glass may reach equilibrium after many 
thousands of years, it will be glassy on the timescale of a 
few years. In such non-equilibrium systems the elapsed 
time t , sometimes called the observation time, waiting 
time, or age, is an important variable for the descrip- 
tion of the macroscopic state of the system. In realistic 
glass forming liquids several types of slow processes may 
occur. Here we shall restrict ourselves to model glasses 
with one type of slow process, because this already ex- 
plains the basic problem. For such systems we have found 
that time shows up through one or two time-dependent 
parameters, namely the effective temperature Te{t) 
and possibly the effective pressure Pe (t) ||lj] . 

Te was introduced phenomenologically half a century 
ago by Tool |3l|] to model the shapes of the specific heat 
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curves on cooling and heating through the glass transi- 
tion. This effective or fitctive temperature describes the 
relaxation of slow modes. As in a plasma, slow and fast 
modes will equilibrate at their own specific temperature. 
We have recently proposed a set of non-linear cooling 
and heating experiments, that should lead to system- 
independent shapes of the specific heat near the glass 
transition ||l|l 

The notion of an effective pressure in the glassy state 
was put forward by Goldstein and Jackie |Q . They 
considered glass formation under pressure. If pressure is 
partly released after the glass has formed, the glassy state 
will depend not so strong on the new pressure, but more 
on its formation pressure, which thus acts as a Active 
pressure. In an exactly solvable model glass we found a 
range of parameters where such a parameter has to be in- 
troduced for a proper thermodynamic description jr^ , p^ . 
The approach with two effective variables is promising, 
since for glass forming liquids it was found experimen- 
tally that at given T,p the same volume can be reached 
via different histories. As these histories then lead to dif- 
ferent futures, the glassy state of these systems cannot be 
coded using only one extra variable Here we shall, 
however, avoid the technical complications of having to 
deal with two effective parameters. 

Let us explain how the effective temperature arises in 
our approach. Consider an aging experiment at fixed 
(T, p) at some long observation time t. Fast modes, with 
relaxation time much less than t, will be in equilibrium 
with the bath temperature T. Modes with timescale 
much larger than t are just quenched. The active slow 
modes are those that have a time scale comparable with 
t. For dynamics on this timescale one may sum out 
the fast processes, by calculating the restricted partition 
sum, i.e. the sum over the fast modes of the Boltzmann 
weights exp(— Ti/fcT), where Ti is the Hamiltonian, i.e. 
the energy of the configuration, and k Boltzmann's con- 
stant. We have worked this out in detail for a model sys- 
tem , but here we donot restrict ourselves to a spe- 
cific system. (How to divide between fast and slow modes 
in realistic systems is a subject of current research). The 
partition sum leads to an effective Hamiltonian of the ac- 
tive slow modes, HoS (active) = U — TScp where U is the 
energy, averaged over the equilibrium processes, and Sep 
their entropy. U and 5'ep depend explicitly on the actual 
configuration of active slow modes, but also on T, which 
results from the fast modes. The effective Hamiltonian 
governs the active slow dynamics, yielding the dynam- 
ical solution for U{t). One may describe this result by 
quasi-static expressions that follow from a partition sum 
at an effective temperature Tg. It is introduced by now 
taking as Boltzmann factor: exp(— Ticff /fcTe). The sum 
over slow modes will yield a generalized free enthalpy 

G{T,T„p)^U~TS,p-T,I + pV (1) 

now U and 5'ep have been averaged over the active slow 
modes, and X is their configurational entropy. It is also 



called information entropy or complexity. T(,(t;T,p) is 
now fixed by matching U{T,Te,p) with the dynamical 
value U{t;T,p). 

In an aging experiment at fixed (T,p) both Tg will 
be a function of time. In a cooling experiment at fixed 
pressure, T and will be functions of time, so that a 
line in the {T,Te,p) space is singled out. For a set of 
smoothly related cooling experiments (e.g. having com- 
mon cooling rates) at different pressures, this leads to 
thermodynamics confined to a surface Te(T, j>) in (T, Tg, 
p) space Q. To cover full space many sets of cooling 
experiments would be needed, e.g. at different cooling 
rates. The results should coincide with other measure- 
ments, such as heating, compression and aging experi- 
ments. 

Very recently numerical data of the (short-time) vibra- 
tional properties in the glassy phase of a binary Lennard- 
Jones system were observed to be governed by an effective 
temperature, in full harmony with our picture Is^j. 



The first Ehrenfest relation 

The two-temperature approach explains in one stroke 
the confusion about the first Ehrenfest relation, eq. (||), 
believed to be never satisfied. The misunderstanding 
arose because in the glassy state the compressibility is 
not a unique material parameter: this measure of the sys- 
tem's response to a change in pressure depends on how 
the measurement is done, and basically on how much 
time the system is given. This fact has already severely 
complicated its experimental determination and interpre- 
tation. In the past, one has typically obtained a and 
dpg/dT from the cooling curves at fixed pressures, but 
determined k by an alternative experiment, such as the 
speed of sound, a procedure which normally works satis- 
factorily. In retrospective, this unvalid approach can be 
traced back to a far too large trust in theories which view 
glass transitions as thermodynamic phase transitions |3^ ] 
prt I pst 0] [|6| . In such a description a unique long-time 
limit of the compressibility exists, that should be mea- 
surable in various ways. McKenna has stressed that in 
experiments on glasses the isothermal compressibility dif- 
fers from the isochoral compressibility . Thus alter- 
native experiments are not allowed, and there remains 
only a tautology to verify. The physical reason is that 
the compressibility strongly depends on the amount of 
time (and other details of the history) that the system 
has to relax. This effect immediately sets in below the 
glass transition. In our approach the validity of the first 
Ehrenfest relation can be checked explicitly 

The second law and the second Ehrenfest relation 

Thermodynamics amounts to giving universal relations 
for the system's state variables at nearby points in (T, 
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Te, p) space. The property dG = -S'cpdT - JdTe + Vdp 
implies together with eq. (|l]) and the first law dU = 
dQ + SW = dQ — pdV^, that the supplied heat is 



dQ = Td5cp + TedJ 



(2) 



The entropy terms are the same as for a cup of coffee in 
a room: also then there are two temperatures, the one 
of the room and of the cup, and there are two entropies. 
On top of that, the timescale for heat exchange is large, 
allowing maintenance of a spontanous difference in tem- 
peratures. A very similar situation occurs in a glass. 
The fast modes (more precisely: fast modes combined 
with the bath) play the role of the bath, while the slow 
configurational modes play the role of the coffee. 

We may decompose the change AS of the total entropy 
S — Sep +1 as the sum of the externally supplied entropy 
deS —dQ/T and the internally produced entropy 



T — T 
d^S - dJ 



(3) 



It is positive because heat flows from high to low temper- 
atures. (At this point plays the role of a real second 
temperature!) The entropy production per unit of time 
is diS/dt. 

Within a set of smoothly related cooling experiments 
one is confined to a surface Tg = Ti,(T,p), pe = pi,{T,p). 
This implies that G has a discontinuous slope at the glass 
transition. This difference in slope was first emphasized 
for spin glasses |3|]. Since its naive explanation as 
being caused by a latent heat does not apply, since that 
has typically not been observed, we were motivated to 
find the physical background, which turned out to be the 
two-temperature approach, discussed here. The differ- 
ence in slopes implies for the Maxwell relation 



U' + pV + TV ^ (Te - Tfe)I' + T±Tl 



(4) 



where we denoted 9/i9tIp by a dot and d/dp\T by a 
prime. In equilibrium (Te = T) the right hand side van- 
ishes, but out of equilibrium there is no reason why it 
should. The second Ehrenfest relation is also modified, 
since it employs this Maxwell relation. It now becomes 



ACp ^ dp 1 - Te dX 

^ = Aa— + — 

TV dT V dT 



(5) 



where d/dT = d/dT + {dpg/dT)d/dp yields the total 
change along the glass transition line. The last two terms 
are new, but would vanish in equilibrium. The violation 
of the second Ehrenfest relation thus gives information 
about the configurational entropy. 

For atactic polystyrene Rehage and Oels report a value 
n = 1.09 K, 1. Having realized that these authors used a 
short-time value for the compressibility in the glass, we 
have taken the theoretical value from the first Ehrenfest 
relation, and obtained 11 = 0.77 This can be de- 

scribed by eq. (^). Values 11 < 1 also occur in the exact 
dynamical solution of a toy model [O. 



In the past it was proven that 11 > 1 is needed for 
mechanical stability ||27| . The starting point was the as- 
sumption that at the glass transition unspecified order 
parameters undergo a thermodynamic phase transition. 
It was even shown that the same assumptions implied the 
equality 11 = 1 [^6| . Since experiments mostly yielded 
2 < n < 5, this remained a field of confusion. It is now 
clear that the assumption of a thermodynamic transition 
is false: the above theoretical framework does not involve 
it, and the Rehage-Oels experiment has 11 < 1. 

A phase diagram of amorphous ice was presented re- 
cently | |40[ |. This approach, however, employs equilibrium 
thermodynamics. In particular, it assumes the validity of 
the Clausius-Clapeyron relation at the first order tran- 
sition from low density amorphous ice to high density 
amorphous ice. However, for first order transitions to a 
glassy state the Clausius-Clapeyron relation will be mod- 
ified icil 



dp, ^U + pAV dUT,p,{T)) 
^ dT~ T +^ dT > ' 

It would be interesting to test this for amorphous ice 



Apparent violations of the laws of thermodynamics 

The basis of our approach has been to show that glassy 
systems respect the fundamental laws of thermodynam- 
ics. However, if one "forgets" that the material is in 
a non-equilibrium state, glass addresses practically all 
thermodynamic paradoxes that have been discussed in 
the past. 

Take a piece of natural or volcanic glass and keep it at 
room temperature for extremely long time. The material 
is in an excited, non-equilibrium state. In the course of 
time it will slowly evolve towards the equilibrium crys- 
talline state. In doing so it will go to a state of lower 
energy and configurational entropy. As is well known, 
see e.g. and described by eq. (||), this will gener- 

ate heat. There will be no true violation of the first and 
second law: after the glass has reached equilibrium, such 
processes are impossible. 

The released heat can be used in several ways: 



Energy out of "nothing" : glass as a source of energy 

An efficient way to extract all configurational energy 
of a glass is as follows. One first brings it into its liq- 
uid state by heating; then one lets it carefully go into 
the crystalline state. This process allows to extract the 
configurational energy in a reasonable time. The energy 
needed for heating can in principle be regained upon cool- 
ing. On top of that, one obtains the crystallization heat 
AU (also called relaxation heat or latent heat), corrected 
for the difference in specific heats of glass and crystal. 
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For metallic glasses this principle is well known to ma- 
terial scientist |4^, and it equally holds for glasses and 
other amorphous materials, such as glassy volcanic rocks. 
In analogy with metal alloys we can estimate the con- 
figurational energy as some 10 kJ /mole atoms. This is 
some 30 times less than the conventional gain of energy 
from combustion of oil and stone coal. This principle 
may also be applied to extract the energy due to stresses 
in polycrystalline rocks. Though these mechanisms are 
presently not attractive, they may be relevant for travel- 
ling or living in outer space, or perhaps for the mining of 
asteroids 1431. 



to zero, in particular in cases where the glass formation 
temperature is much larger than the room temperature. 

It should be admitted that for a small piece of glass 
such a mobile is merely an academic object. In the above 
terrestrial analog, the heat released by radioactive decay 
finally performs work by moving of the earth crust. 



The third law and the residual entropy 

It is long known that on realistic time scales the zero 
temperature entropy of glass does not vanish. It is equal 
to the confiurational entropy of slow modes I. 



Heat out of "nothing": an isolated piece of glass heats itself 
slowly up 



The zeroeth law and the effective temperature 



When the glass is kept in an isolated container (dQ = 
0) , the heat will be absorbed by the fast processes of the 
glass. As they determine the temperature, the glass will 
slowly heat up. Eq. (||) says that TScp + Tj = 0. In the 
equilibration process will go down, but T, the tem- 
perature that can be read off from a thermometer, will 
enhance. 

An analogous manifestation is the fact that most of the 
earth's interior is fluid due to radioactive decay. Indeed, 
the radioactive content, with its decay time of several bil- 
lion years, is a slow and thus "configurational" process. 
The generated heat is absorbed by the "fast" geological 
processes, which maintain the fluid nature. 



The configurational modes of a window glass "live" at 
an effective temperature. During the glass formation pro- 
cess they got essentially stuck after the glass transition, 
which for window glass occurs at some 1000 K . As they 
did not relax very much since then, the effective tem- 
perature Te of these modes is still close to, but below, 
the glass transition temperature. It is thus not proper to 
say that a glass has a single, constant temperature: the 
one related to its slow modes is very different from room 
temperature, and will be measured when performing ex- 
periments on long enough timescales. 



APPENDIX: THE EHRENFEST RELATIONS 



The guantum glass as an almost ideal motor 

A quantum glass in a bath at very low temperature, 
or a realistic glass in outer space with its 3 Kelvin back- 
ground temperature, can have a high effective tempera- 
ture, but basically have vanishing bath temperature. The 
efficiency for transforming the dynamically released heat 
into work, rj = 1 — T/T,,, can be arbitrarily close to unity, 
at the verge of violating the second law. 



Work out of "nothing": a "perpetuum mobile of the third 
kind" 



The heat released by the glass can also be used to per- 
form work by a suitable machine. We thus have a mobile 
that performs work as long as the glass is not in equilib- 
rium. From geology it is known that, at least for Si-rich 
glasses, this period may be astronomical [|l|. We may 
therefore call the non-equilibrium glass, combined with 
some apparatus that performs the work, a "perpetuum 
mobile of the third kind". Indeed, work dW ~ — dQ is 
performed at the cost of configurational entropy, which, 
according to the third law, should have been much closer 



At a first order phase transition from phase A to 
phase B, the slope of the transition line is given by the 
Clausius-Clapeyron relation. It derives from continuity 
of the free enthalpy G = U + pV - TS. The difference 
AG = Ga — Gb vanishes along the transition line. Dif- 
ferentiating AG{T,pg{T)) — one obtains 



AU + pAV ^ dpg 
T dT 



(7) 



For second order phase transitions of classical type one 
can differentiate AV = 0, which yields the first Ehrenfest 
relation 



Aa = An 



dT 



(8) 



It is said to be violated in glasses, but we shall argue 
against that. 

By considering AU = one obtains the second Ehren- 
fest relation 



^=Aa^ 
TV dT 



(9) 



It is believed to be often satisfied for glasses, though not 
in all cases. We shall show that in principle it is violated. 

The information of both Ehrenfest relations is often 
coded in the so-called Prigogine-Defay ratio 
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TV{Aa)^ ^ ' 

Using eq. (H) and (^) it seen to equal to unity in equilib- 
rium. For glasses it is reported to lie typically between 
2 and 5. We shall stress that also values below unity are 
possible, and actually already have been measured. 
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FIG. 1. Data of the glass transition for cooling atactic polystyrene with 18 K/h, taken from Rchagc and Ocls (1976): specific 
volume V [cm^ /g) versus temperature T (K) at various pressures p (kbar). The data in the liquid essentially lie on a smooth 
surface, and so do the data in the glass. The first Ehrenfest relation describes no more than the intersection of these surfaces, 
and is thus automatically satisfied. The values for the compressibility that follow from these data will generally differ from 
results obtained via other experiments. 
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